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Abstract 

In a recent formulation of a quantum field theory of forward rates, 
the volatility of the forward rates was taken to be deterministic. The 
field theory of the forward rates is generalized to the case of stochastic 
volatility. Two cases are analyzed, firstly when volatility is taken to be 
a function of the forward rates, and secondly when volatility is taken 
to be an independent quantum field. Since volatiltiy is a positive 
valued quantum field, the full theory turns out to be an interacting 
nonlinear quantum field theory in two dimensions. The state space 
and Hamiltonian for the interacting theory are obtained, and shown to 
have a nontrivial structure due to the manifold moving with a constant 
velocity. The no arbitrage condition is reformulated in terms of the 
Hamiltonian of the system, and then exactly solved for the nonlinear 
interacting case. 

PACS:02.50.-r Probability theory, stochastic processes 05.40.+j 
Fluctuation phenomena, random processes and Brownian motion 03.05.- 
w Quantum mechanics 
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1 Introduction 



Forward rates are essential to the debt market, and have wide-ranging ap- 
phcations in finance. The most widely used model of the forward rates is 
the the Heath- Jarrow- Morton (HJM) There are a number of ways that 
the HJM-model can be generalized. In references P| and P] a correlation 
between forward rates with different maturities was introduced, and in 0] 
and the forward rate was modeled as a stochastic string. 

The application of techniques of physics to finance 0, ^ have proved 
to be a fruitful field; in particular path integral techniques [|| have been 
applied to various problems in finance. In P] path integral techinques were 



applied to study the case of a security with stochastic volatility. In |10 the 



HJM-model was generalized by treating the forward rates as a quantum field; 
empirical studies in |Tl| show that the field theoretic model for the forward 
rates proposed in [l^ fits market data fairly well. 

The volatility of the forward rates is a central measure of the degree to 
which the forward rates fiuctuate. In the model studied in [1^, the volatility 



of the forward rates was taken to be deterministic. The question naturally 
arises as to whether the volatility itself should considered to be a randomly 
fiuctuating quantity. The volatility of volatility is an accurate measure of the 
degree to which volatility is random. Market data for the Eurodollar futures 
provides a fairly accurate estimate of the forward rates for the US dollar, and 
also yields its volatility of volatility of the forward rates. 

Eurodollar futures data given in Figure |l| plots the 30 days moving av- 
erage of the volatility of volatility for the forward rates, and shows that it 
contributes about 0.0006 — 0.0007 per year to the forward rates, which are in 
the neighbourhood of 0.05 — 0.06 per year. The fiuctuations in the volatility 
of the forward rates are about 10% of the forward rates, and hence significant. 

We conclude from the data that the volatility of the forward rates needs 
to be treated as a fiuctuating quantum (stochastic) field. The widely stud- 
ied HJM-model has been further developed by |]T^ (and references cited 
therein) to account for stochastic volatility. Amin and Ng [jl3| studied the 
market data of Eurodollar options to obtain the implied forward rate volatil- 



ity and Bouchaud et. al. [14] analyzed the future contracts for the forward 



rates. Both references concluded that many features of the market, and in 
particular the (stochastic) volatility of forward rate curve, could not be fully 
explained in the HJM-framework. 

The model for the forward rates proposed in is a field theoretic gen- 
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Figure 1: The Volatility of the Volatility of the Forward Rates 



eralization of the HJM-model, and so it is natural to extend the field theory 
model to account for stochastic volatility of the forward rates. In contrast 
to quantum field theory, the formulation of the forward rates as a stochastic 
string cannot be extended to the case where volatility is stochastic due 

to nonlinearites inherent in the problem. 

The forward rates are the collection of interest rates for a contract entered 
into at time t for an overnight loan at time x > t. At any instant t, there 
exists in the market forward rates for a duration of TpR in the future; for 
example, if t refers to present time to, then one has forward rates from to till 
time to + TpR in the future. In the market, TpR is about 30 years, and hence 
we have Tp^ > 30 years. In general, at any time t, all the forward rates exist 



till time t + Tpn [IC]. The forward rates at time t are denoted by f{t,x), 
with t < X < t + TpK, and constitute the forward rate curve. 

Since at any instant t there are infinitely many forward rates, it resembles 
a (non-relativistic) quantum string. Hence we need an infinite number of 
independent variables to describe its random evolution. The generic quantity 
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describing such a system is a quantum field [T^. For modeling the forward 



rates and Treasury Bonds, we consequently need to study a two-dimensional 
quantum field on a finite Euclidean domain. 

We consider the forward rates f{t, x) to be a quantum field; that is, 
f{t, x) is taken to be an independent random variable for each x and 
each t. For notational simplicity we consider both t and x continuous, and 
discretize these parameters only when we need to discuss the time evolution 
of the system is some detail. 

In Section 2, we briefly review the quantum field theoretic formulation 
given in [|I^ of the forward rates with deterministic volatilty. In Section 3 



the case of stochastic volatility is anlayzed, and which can be done in two 
different ways. Firstly, volatility can be considered to be a function of the 
(stochastic) forward rates, and secondly volatility can be considered to be 
an independent quantum field. Both these cases are analyzed. The resulting 
theories are seen to be highly non-trivial non-linear quantum field theories. 

In Section 4 the underlying state space and operators of the forward 
rates quantum field is defined. In particular the generator of infinitesimal 
time evolution of the forward rates, namely the Hamiltonian, is derived for 
the two cases of stochastic volatility. In Section 5 the Hamiltonian for the 
forward rates with stochastic volatility is derived. In Section 6 a Hamiltonian 
formulation of the condition of no arbitrage is derived. In Section 7 the no 
arbitrage constraint for the case of stochastic volatility is solved exactly using 
the Hamiltonian formulation. And lastly, in Section 8 the results obtained 
are discussed, and some remaining issues are addressed. 



2 Lagrangian for Forward Rates with Deter- 
ministic Volatility 

We first briefly recapitulate the sahent features of the field theory of the 
forward rates with deterministic volatility |Tp . 

For the sake of concreteness, consider the forward rates starting from 
some initial time Tj to a future time t = Tf. Since all the forward rates 
f{t, x) are always for the future, we have x > t; hence the quantum field 
f{t,x) is defined on the domain in the shape of a parallelogram V that is 
bounded by parallel lines x = t and x = TpR + 1 in the maturity direction , 
and by the lines t = Ti and t = Tf in the time direction, as shown in Figure 
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(|^). Every point inside the domain V represents an independent integration 
variable /(t, x). 

The field theory interpretation of the evolution of the forward rates, as 
expressed in the domain P, is that of a (non-relativistic) quantum string 
moving with unit velocity in the x (maturity) direction. 

A 

t 




Figure 2: Domain V of the Forward Rates 



Since we know from the HJM-model that the forward rates have a drift 
velocity a{t,x) and volatitility these have to appear directly in the 

Lagrangian for the forward rates. To define a Lagrangian, we firstly need a 
kinetic term, denoted by /^kinetic, that is necessary to have a standard time 
evolution for the forward rates. 

We need to introduce another term to constrain the change of shape of 
the forward rates in the maturity direction. The analogy of this in the case 
of an ordinary string is a potential term in the Lagrangian which attenuates 
sharp changes in the shape of the string, since the shape of the string stores 
potential energy. To model a similar property for the forward rates we cannot 

use a simple tension-like term (7^)^ in the Lagrangian since, as we will show 

ox 

in Section 7, this term is ruled out by the condition of no arbitrage. 

The no arbitrage condition requires that the forward rates Lagrangian 

contain higher order derivative terms, essentially a term of the form ( o 

(y Jb V 

such string systems have been studied in |]16| and are said to be strings with 
finite rigidity. Such a term yields a term in the forward rates Lagrangian, 
namely frigidity, with a new parameter fj,; the rigidigy of the forward rates 
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is then given by — and quantifies the strength of the fluctuations of the 

forward rates in the time-to-maturity direction x. In the limit of — > 0, we 
recover (upto some rescahngs) the HJM-model, and which corresponds to an 
infinitely rigid string. The action for the forward rates is given by 

S[f] = dt dxC[f] (1) 

Jt, Jt 

m (2) 

V 

with the Lagrangian density C[f] given by 

^[f] = -^^kineticf/] + -^^rigidityl/] (3) 

J u?\dx\ 



2 It a{t,x) J ^Adx\ a{t,x) )) . 
-oo < f{t,x) < +00 (4) 

The presence of the second term in the action given in eq.(|^) is not ruled 
out by no arbitrage , and an empirical study provides strong evidence 



for this term in the evolution of the forward rates. 

In summary, we see that the forward rates behave like a quantum string, 
with a time and space dependent drift velocity a{t, x), an effective mass given 

by — -, and string rigidity proportional to — . 

a{t,x) jji^ 

Since the field theory is defined on a finite domain V as shown in Figure. 
H, we need to specify the boundary conditions on all the four boundaries of 
the finite parallelogram V. 

• Fixed (Dirichlet) Initial and Final Conditions 

The initial and final (Dirichlet) conditions in the time direction are 
given by 

T,{Tf)< X <T,{Tf)+TFR : f{T„x), f{Tf,x) (5) 
: specified initial and final forward rate curves 

• Free (Neumann) Bounday Conditions 

To specify the boundary condition in the maturity direction, one needs 
to analyze the action given in eq.(|lD and impose the condition that 
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there be no surface terms in the action. A straightforward analysis 
yields the following version of the Neumann condition 



tS \it.J) ) = ° 

: x = toTX = t + TpR (7) 

The quantum field theory of the forward rates is defined by the Feynman 
path integral by integrating over all configurations, and yields of /(t, x) and 

Z = ! Dfe^^^^ (8) 



Df ^ H dfit^x) (9) 

{t,x)eV 

Note that c^t-^'/ Z is the probability for different field configurations to 
occur when the functional integral over f{t, x) is performed. 



3 Lagrangian for Forward Rates with Stochas- 
tic Volatihty 

To render the volatility function a{t, x) stochastic, in the formalism of quan- 
tum field theory, requires that we elevate a{t, x) from a deterministic function 
into random function, namely into a quantum field. There are essentially two 
ways of elevating volatility to a stochastic quantity, namely to either (a) con- 
sider it be a function of the forward rate f(t, x), or else (b) to consider it to be 
another independent quantum field a{t, x). We study both these possibilities. 



3.1 Volatility a Function of the Forward Rates 

We consider the first case where volatility is rendered stochastic by making it 
a function of the forward rates [|l3l . The standard models using this approach 



consider that volatility is given by 



ait, X, fit, x)) = aoit,x)rit,x) (10) 
cro(t,x) : Deterministic function (11) 
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Since volatility a{t,x) > 0, we must have f{t,x) > 0. Hence, in contrast 
to eq.(^), we have 

f{t, x) = /oe'^(*'^) > ; - oo < x) < +oo (12) 

Having f{t, x) > is a major advantage of the model since in the financial 
markets forward rates are always positive. In the limit of // — > 0, the following 
HJM- models are covered by eq.(|l^), and these models have been discussed 



from an empirical point of view in |13 



1. Ho and Lee (1986) Model : a{t, x, f{t, x)) = ao 

2. CIR (1985): a{t, x, f{t, x)) = aof^t, x) 

3. Courtadon(1982) -.aitjX, f{t,x)) = aof{t,x) 

4. Vasicek (1977) :cr(t, x, /(t, x)) = ctq exp(— A(x — t)) 

5. Linear Proportional HJM (1992) : a{t,xj{t,x)) = [ao+ai{x-t)]f{t, x)] 

How do we generalize the Lagrangian given in eq.(^ to case where the 
forward rates are always positive? We interpret the Lagrangian given in 
eq.(^ to be an approximate one that valid only if all the forward rates are 
close to some fixed value /q. We then have 



at /«^*""^ (13) 

/o^ + 0(^') (14) 



Hence we make the following mapping 



dt dt 

Eq.(^ then generalizes to 



df{t,x) ,d(p{t,x) , 
' Jo — TT, — 



C[(p] = /^kinetic M + -^^rigidityl^] (16) 

■^mA -a(t,x) Y ^ 1 J- d / /o^%^ - Q(t,a;) \ 1 
cro(t,x)e^*(*'^) J n'^XdxK ao(t, x)e'^<^(*'^) // 
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We will show later - in deriving the Hamiltonian - that the system needs a 
non-trivial integration measure. We hence define the theory by the Feynman 
path integral 



f^^sm (17) 

/r+oo 
Dcj^r = n / dct^it.x) r{t,x) (18) 



{t,x)eV ■ 



The boundary conditions given for f{t, x) in eqs.(|^) and (|^) continue to hold 
for stochastic volatility Lagrangian given in eq. (|l6|) . 



3.2 Volatility as an Independent Quantum Field 

We consider the second case where volatility a{t, x) is taken to be an inde- 
pendent (stochastic) quantum field. Since one can only measure the effects 
of volatility on the forward rates, all the effects of stochastic volatility will 
be manifested only via the behaviour of the forward rates. 

For simplicity, we consider the forward rate to be a quantum field as given 
in eq.(|D with 

f{t,x): - OO < /(t, x) < +00 (19) 

Since the volatility function a{t, x) is always positive, that is, a{t, x) > 0, 
we introduce an another quantum field h{t, x) by the following relation (the 
minus sign is taken for notational convinience). 

a{t,x) = aoe''''^^'''\ - oo < h{t, x) < +oo (20) 

The system now consists of two interacting quantum fields, namely 
f{t,x) and h{t,x). The interacting system's Lagrangian should have the 
following features. 

• A parameter ^ that quantifies the extent to which the field h{t, x) is 
not deterministic. A limit of ^ would, in effect, 'freeze' all the 
fluctuations of the fleld h{t, x), and reduce it to a deterministic function. 

• A parameter k to control the fluctuations of h{t, x) in the maturity 
direction similar to the parameter /i that controls the fluctuations of 
the forward rates f{t, x) in the maturity direction x. 
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• A parameter p with — 1 < p < +1 that quantifies the correlation of the 
forward rates' quantum field f{t,x) with the volatility quantum field 
h{t, x). 

• A drift term for volatility, namely l3{t,x) - which is analogous to the 
drift term a{t,x) for the forward rates. 

The Lagrangian for the interacting system is not unique; there is a wide 
variety of choices that one can make to fuUfil all the conditions given above. 
A possible Lagrangian for the interacting system, written by analogy with 
the Lagrangian for the case of stochastic volatility for a single security |^ , is 
given by 



2(1- 






' d 


2/.2 ( 


dx 




(21) 



S[f, h]= C (22) 



with action 



r 

We need to specify the boundary conditions for the interacting system. 
The initial and final conditions for the forward rates f{t,x) given in eq.(^ 
continue to hold for the interacting case, and are similarly given for the 
volatility field as the following. 

• Fixed (Dirichlet)Initial and Final Conditions 

The initial value is specified from data, that is 

Ti{Tf)< X <Ti{Tf)+TFR, a{Ti,x),a{Tf,x) (23) 
: specified initial and final volatility curves 



The boundary condition in the x— direction for the forward rates f{t, x) as 
given in eq.(^ - continues to hold for the interacting case, and for volatility 
field is similarly given by 



10 



• Free (Neumann) Boundary Conditions 

: X = t or X = t + Tpji (25) 

On quantizing the volatility field a{t,x) the boundary condition for the for- 
ward rate f{t, x) given in eq.(|^) is rather unusual. On solving the no arbitrage 
condition, we will find that a is a (quadratic) functional of the volatility field 
a{t, x); hence the boundary condition eq.(^ is a form of interaction between 
the f{t, x) and a{t, x) fields. 

We need to define the integration measure for the quantum field h{t,x); 
the derivation of the Hamiltonian for the system dictates the following choice 
for the measure, namely 



/n+oo 
DfDa-^ = n / df{t,x)da~\t,x) (26) 

/ + 00 
df{t,x)dh{t,x)e''^''^^ (27) 
... -oo 



(t,x)eP 



The partition function of the quantum field theory for the forward rates 
with stochastic volatility is defined by Feynman path integral as 

Z= [ DfDa-^e^^f^^^ (28) 



The (observed) market value of a financial instrument, say /i], is 
expressed as the average value of the instrument - denoted by < 0[f, h\ > 
- taken over all possible values of the quantum fields /(t, x) and h{t, x), with 
the probability density given by the (appropriately normalized) action. In 
symbols 

< 0[f, h]>=^l DfDa-' 0[f, /i]e^[^''^l (29) 

We consider the limit of the volatility being reduced to a determinis- 
tic function. For this limit we have ^, p and k — > 0. The kinetic term of 
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the h{t,x) field in the action given in eq.(P^ has the hmit (upto irrelevant 
constants) 

>'5„n-p{4X(^)Vn^(f-'') (30) 

which implies that 

< a{t, x)> = (To < e-'^(*'^') > (31) 
= aoexp{- [ dt'P{t',x)} + 0{C,K,p) (32) 

J to 

4 Hamiltonian and State Space 

The Feynman path integral formulation given in eqs. (|T^ and ( pSj ) is useful 
for calculating the expectation values of quantum fields. To study questions 
related to the time evolution of quantities of interest, one needs to derive the 
Hamiltonian for the system from its Lagrangian. Note the route that we are 
following is opposite to the one taken in ^ where the Lagrangian for a stock 
price with stochastic volatility was derived starting from its Hamiltonian. 
The state space of a field theory is a hnear vector space - denoted by V 

- that consists of functionals of the field configurations at some fixed time t. 
(A brief discussion of the state space is given in |^.) The dual space of V 

- denoted by Vouai ^ consists of all linear mappings from V to the complex 
numbers, and is also a linear vector space. Let an element of V be denoted 
by \g > and an element of Vouai by < p|; then < p\g > is a complex number. 
We will refer to both V and Vouai as the state space of the system. The 
Hamiltonian Ti. is an operator - the quantum analog of energy - that is an 
element of the tensor product space V ® VbubI- The matrix elements of H 
are complex numbers, and given hj < p\H\g >. 

In this section, we study the features of the state space and Hamiltonian 
for the forward rates. For notational brevity, we consider the forward rates 
quantum field f{t, x) to stand for both the quantum fields f{t, x) and h{t, x). 
Since the Lagrangian for the forward rates given in eq. (pTD has only first order 
derivatives in time, an infinitesimal generator, namely the Hamiltonian Ti 
exists for it. Obtaining the Hamiltonian for the forward rates is a complicated 
exercise due to the non-trivial structure of the underlying domain V. In 
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particular, the forward rates quantum field will be seen to have a distinct 
state space Vj for every instant t. 

For greater clarity, we discretize both time and maurity time into a finite 
lattice, with lattice spacing in both directions taken to be e. (For a string 
moving with velocity v, the maturity lattice would have spacing of ve.) On 
the lattice, the minimum time for futures contract is time e; for most appli- 
cations e =one day. The points comprising the discrete domain V are shown 
in Figure |^. 



Figure 3: Lattice in Time and Maturity Directions 



The discrete domain V is given by 

(t, x) e(n, /) ; n,l : integers (33) 

{Ti,Tf,Tpn) e{N,,Nf,NpR) (34) 

Lattice V = {{n, l)\Ni < n < Nf ; n < I < {n + Nfr)} (35) 

f{t,x) ^ fn,i (36) 

df{t,x) ^ fn+l,l - fn,l _ df{t,x) _ fn,l+l - fn,l 

dt e ' dx e 

The partition function is now given by a finite multiple integral, namely 



Z 
S 



{n,l)eV 



S[f\ 



(38) 
(39) 



Consider two adjacent time slices labelled by n and n + 1, as shown in 
Figure ^. S{n) is the action connecting the forward rates of these two time 
slices. 
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Figure 4: Two Consecutive Time Slices for t = ne and t = [n + l)e 



As can be seen from Figure for the two time slices there is a mismatch 
of the 2-lattice sites on the edges, namely, lattice sites (n, n) at time n and 
(n + l,n + 1 + Nfr) at time n + 1 are not in common. We isloate the 
un-matched variables and have the following 

Variables at time n : 

{/„,„, Fn} ; Fn = {fn,i\n <l<n + Nfr} (40) 
Variables at time (n + 1) : 

{Fn, fn+l,n+l+NFR 

Note that although the variables F„ refer to time n + 1, we label it with 
earlier time n for later convinience. From Figure ^ we see that both sets of 
variables F^ and F„ cover the same lattice sites in the maturity direction, 
namely n + 1 < / < n + Nfr, and hence have the same number of forward 
rates, namely Nfr — 1- The Hamiltonian will be expressable solely in terms 
of these variables. 

From the discretized time derivatives defined in eq.(0) the discretized 
action S{n) contains terms that couple only the common points in the lattice 
for the two time slices, namely the variables belonging to the sets F^, Fn- We 
hence have for the action 

S{n) = e 5^ /„,+i,i] (42) 

{1} 

= eJ2^n[Fn;Fn] (43) 
{0 
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As shown is in Figure the action for the entire domain V shown in 
Figure |] can be constructed by repeating the construction given in Figure | 
and summing over the action S{n) over all time Ni < n < Nf. 




Figure 5: Reconstructing the Lattice from the Two Time Slices 

The Hamiltonian of the forward rates is an operator that acts on the 
state space of states of the forward rates; we hence need to determine the 
co-ordinates of its state space. 

Consider again the two consecutive time slices n and n+1 given in Figure 
|[ We interpret the forward rates for two adjacent instants, namely {/n,n, F^} 
and {Fn, fn+i,n+i+NpR} given in eq . (|40|) - and which appear in the action 
eg. (^2]) - as the co-ordinates of the state spaces Vn and V„+i respectively. 

For every instant of time n there is a distinct state space V„, and its dual 
Vouai.n- The co-ordinates of the state spaces V„ and V„+i are given by the 
tensor product of the space of state for every maturity point /, namely 

<fn\ = (g) < fnA =< fn,n\ < Fn\ (44) 

n<l<n+Npj{ 

: CO — ordinate state of VDuai.n 

l/n+l > = \fn+l,l >=\Fn> \ fn+l,n+l+Npn > (45) 

{n+l)<l<n+l+Npii 

: CO — ordinate state of Vn+i 

The state vector |F„ > belongs to the space space V„+i, but we reinterpret 
it as corresponding to the state space JF„ at earlier time n. This interpre- 
tation allows us to study the system instantaneously using the Hamiltonian 
formalism. 
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The state space Vn consists of all possible functions of Npji forward rates 
{fn,n, -Pn}- The state spaces Vn differ for different n by the fact that a different 
set of forward rates comprise its set of independent variables. 

Although the state spaces V„ and Vn+i are not the identical, there is an in- 
tersection of these two spaces, namely V„nV„+i that covers the same interval 
in the maturity direction, and is coupled by the action S{n). The intersecion 
yields a state space, namely on which the Hamiltonian evolution of the 
forward rates takes place. In symbols, we have 



Vn+l = ® \fn+l,n+l+NFji > 

"l^Dual.Ti = < fn,n\ ® -^Dual.n 



(46) 
(47) 
(48) 



The Hamiltonian 7i„ is an element on the tensor product space spanned by 
the operators |F„ >< namely the space of operators given by JF„ ® 

•^Dual,rf 

The vector spaces V„ and the Hamiltonian 7i„ acting on these spaces is 
shown in Figure 























\ 





+2 



Vn+l 

v„ 



Figure 6: Hamiltonians Tin propagating the space of Forward Rates V, 



Note that both the states |F„ > and < F„| belong to the same state 
space Tni and we use twiddle to indicate that the two states are different; 
in contrast, for example, the two states |/ > and < /| indicate that one state 
is the dual of the other. 

As one scans through all possible values for the forward rates / and /, one 
obtains a complete basis for the state space Vn- In particular, the resolution 
of the identity operator for V„ - denoted by X„ - is a reflection that the basis 
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states are complete, and is given by 

n<l<n+Npii 

— j dfji^n dFn \fn,m fn,m -^n| (^0) 

The Hamiltonian of the system Ti is defined by the Feynman formula 
(upto a normalization), from eg. (^2]) , as 

^ p'^'E{l}^n[f„,l,fn+l,l] = / f P \P~'''^"\F f ,1 , 1 , l\r > f5ll 

yn'^ — ^ Jn,ni ^ 1-^ ni J n+l,n+l+Npii ^ \^^) 

where in general pn is a field-dependent measure term. Using the property 
of the discrete action given in eq. (|^) , we have 

= < Fn\e-''^"\Fn> (53) 
Equation ( |5^ ) is the main result of this Section. 

In going from eq.([52|) to eq. (^3]) we have used the fact that the action 
connecting time slices n and n + 1 does not contain the variables and 
fn+i^n+i+NpR respectively. This leads to the result that the Hamiltonian Hn 
consequently does not depend on these variables. 

The interpretation of eq.(|53|) is that the Hamiltonian 7i„ propagates the 
initial state < Fn\ in time e to the final state \Fn >■ Note the relation 

>=< F„|e~^^"|F„ > (54) 

shows that there is an asymmetry in the time direction, with the Hamiltonian 
being independent of the earliest forward rate fn^n of the initial state and of 
the latest forward rate fn+i,n+i+NpR of the final state. It is this asymmetry 
in the propagation of the forward rates which yields the parallelogram domain 
V given in Figure 0, and refiects the asymmetry that the forward rates f{t, x) 
exist only for x > t. 

For notational simplicity, we henceforth use continuum notation; in par- 
ticular, the state space is labelled by Vt, and state vector by \ft >. The 
elements of the state space of the forward rates Vj includes all the financial 
instruments that are traded in the market at time t. In continuum notation. 
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from eq.(^), we have that 



\ft> = 09 \f{t,x)> (55) 

t<x<t+Tpji 

\Ft> = (g) \f{t,x)> (56) 

t<x<t+Tpji 

In continuum notation, the only difference between state vectors \ft > and 
I Ft > is that, in eg. (|5^) , the point x = t is excluded in the continuous tensor 
product. 

The partition function Z given in eq.(38) can be reconstructed from the 
Hamiltonian by recursively applying the procedure discussed for the two time 
slices. We then have, in continuum notation, that 



D/e^[^] (57) 

= < /initialirj exp(- f ' n{t) rft)}|/final > (58) 

where the symbol T in the equation above stands for time ordering the (non- 
commuting) operators in the argument, with the earliest time being placed 
to the left. 



4.1 Bond State Vectors 

The most important state vectors in finance are those of the coupon and zero 
coupon bonds. 

Consider a risk-free zero coupon Treasury bond that matures at time T 
with a payoff of $1. The price of the bond at time t < T is given by 

P{t, T) = e- = t, T] (59) 

The state vector |P(t, T) > is an element of the state space V*. We write 
the bond state vector is as follows 

F(t,T) = <ft\P{t,T)> (60) 
= e-5tSi^y^ (61) 

Another state vector is the coupon bond |S >, with payoff of amount q 
at time T/, with a final payoff of L at time T. We then have that the state 
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vector of the coupon bond is linear superpostion of the zero coupon bonds 
and is given by 



\B{t) >=J2ci\Pit,Ti) > +L\P{t,T) > (62) 

5 Hamiltonian for the Forward Rates with 
Stochastic Volatility 

We have obtained a general expression for the Hamiltonian in terms of the 
action S as given in eq.(|^), and need to apply this formula to the case of 
the specific Lagrangian of the forward rates to obtain the explicit expression 
for its Hamiltonian Ti. 

From eg . (1531) we have the following 

p„e^(") = p^e^Eo}^n[^^n,^„;KA] (g3) 
= <F„;^Je-^«"|F„;//„> (64) 

where we have explicitly included the volatility quantum field h{t, x) in the 
equation above. 

For notational simplicity, we consider the maturity direction x to be con- 
tinuous, and consider only the time direction to be discrete. In the continuum 
notation, the subtleties of the variables at time t and t + e are accounted for 
by carefully analyzing the variables appearing on the boundaries of the in- 
terval [t < X < t + Tfr]- We have, for the action S{n) for time t = ne, the 
following 



S{n) = e J Cnit.x) (65) 
dx (66) 

5.1 Hamiltonian for the Forward Rates with Stochas- 
tic VolatiUty 

As a warm-up exercise, we first obtain the Hamiltonian for the simpler case 
of the volatility being a function of the forward rates. Recall the Lagrangian 
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for this system is given by 

C[(j)\ = /^kinetic [0] + -Crigidityl^] (67) 



1 



On discretizing the lagragian we obtain, using boundary condition eq.(|]), 
that 

^(">-/^" = 4/4-^5^)^ (68) 

A = /o<r„"'e-''*(A+, - 0, - e/„-'a) (69) 

We rewrite eq.(p8|) using Gaussian integration and obtain (ignoring hence- 
forth irrelevant constants) 

gS(n) = Y[J dp{x) exp { - ^ y p{x)D{x, x'; t)p{x') + i j p{x)A{x)] (70) 

,2 



1 d' 

[1 -—^)D[x,x'] t) = 5{x — x') ; Neumann boundary conditions 



An exphcit derivation of the propagator D(x, x'\ t) is given in the Appendix, 
and yields 

D{x,x']t) = — -[cosh(//TpR - ^\x-x'\) 
2 smh (/iTp/j) 

+ cosh(//Tpij - fi{x + x' - 2t))] (71) 
Let the measure term be defined as 

Pn = l[r''{ne,x) (72) 

X 

Rescale p{x) as follows (and which will effectively cancell the measure term) 

P /o'Voe'^'^p (73) 

We then have 

p^e-^^"^ = Ylj rfp(x)exp|z j p{x){(l)t+e- <Pt- €fQ^a){x) 



X 

€ 



joMDi.M*^)\ (74) 
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Recall from eq.(p^ that the Hamiltonian is defined by 

p„e^(") = <(t)t\e~''^^\(Pt+e> (75) 



and yields, dropping the subscript t in 0t, the Hamiltonian for the forward 
rates, namely 



(77) 



The Hamiltonian is non-Hermetian with complex eigenvalues. Although this 
would problematic in physics, this is not so in finance since the Hamiltonian 
is not a physical quantity (such as energy) whose eigenvalues are observables, 
and hence is not required to have real eigenvalues. 

5.2 Hamiltonian for the Forward Rates and Volatility 
Quantum Fields 

We now consider the case when both the forward rates and its volatility 
fluctuate independently and are represented by separate quantum fields. We 
hence examine the Lagrangian given in eq.(^T]), namely 

C{t,x) = -— -{^ p^* I dt 



2(1 -p2) I a ^ i 



1 / d (%-A\ 1 / d / t-/9' 



(78) 



where recall 



- oo < f{t,x),h{t,x) < +00 (79) 



Discretizing time, and for notational simplicity suppressing the time and 
maturity labels, we write the Lagrangian C in matrix notation as follows 



-1 J [a-^A ^-^B]ix)Mix,x';t) 



{x') (80) 
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where 



M{x,x';t) 



1 

_ P 



l-p2 

J_8P_ 

1— p2 k2 Q^2 



1 



5{x — x') 



and 



A = ft+e - ft-ea 



(81) 



(82) 
(83) 



Note that in obtaining eq.(pO|) for S{n) we have used the boundary conditions 
on the fields given in eqs.(||) and (^41). 

We rewrite eq.(pOD using Gaussian integration and obtain (ignoring irrel- 
evant constants) 



/ dp{x)dq{x) exp ( — - / [ p q ] A4 ^ 



P 



+i [p q] 



M) 



Recall from eq.(p7D that the Feynman path integral has a non-trivial 
measure cr~^(t,a;), and in obtaining the Hamiltonianwe we need to take this 
into account . 

Define the measure term by 



X 

and rescale the p and q variables in eq. (^4]) for each x as 



p 
q 

We then obtain from eq.(p^) that 

oS(n) 



ap 

eg 



(85) 



(86) 
(87) 



DpDq X 



exp ( - ^ y [ ^q]M ^ 



ap 


-/ 




' f- 


f- 


ea 




[p q] 




h- 


eP _ 



(88) 
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showing that the measure term cancels out. We hence have from above 



5(n) 



<f;h\e-^'^"\f;h> 



(89) 
(90) 



and which yields the Hamiltonian for the forward rates and volatility as 
independent quantum fields given by 



d c 6 
iSf ^ iSh 



M 



-1 



^ iSh 



From eq.(0) we have 

M'^{x, x'; t) = c 
where 



D_-D+ + i^(r+Z}+-r_Z^_ 
p{D- - D+) 



dj oh 



PiD- - D+) 



2(1 - p2) 



(92) 
(93) 



+ f±)D±{x, x'\ t) = 6{x — x) with Neumann boundary conditions 

(94) 



For solving the no aribitrage condition, we will need 

G{x]x',t) = M^l{x,x']t) (95) 

: [k\D_ - D^) + (1 - p^)(r+D+ - r_D_)] 

(96) 



P2 

^(«;2 - n'^y + 4p2p2^2 



6 Hamiltonian Formulation of No Arbitrage 

The principle of no arbitrage is central to the theory of finance, and a path 
integral formulation of this principle is given in [^. For the case of deter- 
ministic volatility, the Lagrangian for the forward rates as given in eq.(^ is 



23 



quadratic, and hence the condition of no abitrage could be solved exactly by 



performing a Gaussian path integration ||10 |. 

For the case of stochastic volatility, the Lagrangian is nonlinear and hence 
the condition of no arbitrage cannot be solved explicitly using the path in- 
tegral; for this reason we reformulate the no arbitrage condition using the 
Hamiltonian. We will show that the Hamiltonian formulation even for the 
nonlinear theory of the forward rates with stochastic volatility allows for an 
exact solution of the no arbitrage condition. 

We first derive the Hamiltonian formulation of no arbitrage for the case of 
a single security S, since the derivation for the forward rates is more complex. 

6.1 No Arbitrage for a Single Security 

Consider an option on a security S = that matures at time T and has a 
payoff function given by g{x, K), where K is the strike price. As discussed in 
[T0| , the risk-free evolution of the security is given by the Hamilitonian Tis, 



with the value of the option at time t < T being given by 

/oo 
dx' < x|e-(^-*)^^|x' >^(x') (97) 
-oo 

(r is a given constant risk-free spot interest rate). 

The martingale condition for the risk-free evolution of the security is that 
the discounted evolution of the future price of the security at some future 
time, say t*, is equal, on the average, to the price of the security at earlier 
time t. The equation for the martingale condition states 

S{x{t)) = E[t,t,][e-^''-'^''S{x{Q)] (98) 

where the notation denotes the average value of Y over all the 

stochastic variables in the time interval (t,t*]. From eq.(^), we have 



/oo 
dx' < x\e'^''-'^'^^\x' > Six') (99) 
-oo 



/oo 
dx' < x|e~(**-*)(^^+'')|x' >< x'\S > (100) 
-oo 

Using the completeness equation for a single security given by 

POO 

1= dx\x >< x\ (101) 
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yields from eq.( |lUU| ), the operator equation 

|5 g-(t*-t)(W.+r)|^ (102) 

Since time is arbitrary, we have 

{ns + r)\S>=0 (103) 

One can easily verify that the Black-Scholes Hamiltonian for both the case 
of deterministic and stochastic volatility given in ||^ satisfies eq. ( |103| ). 

The result given in eq.( p.03| ). shows that the existence of a martingale 
measure is equivalent to a risk-free Hamiltonian that annihilates (upto a 
constant r) the underlying security S. We will see that a similar condition 
holds for the Hamiltonian of the forward rates, but with a number of com- 
plications arising from the non-trivial domain of the forward rates and that 
the spot rate r{t) = f{t,t) is itself a stochastic quantity. 



6.2 No Arbitrage for the Forward Rates 

The principle of no aribitrage states that the price of the bond P{t^,T) at 
some future time T > if:,, > t is the equal to the price of the bond at time t, 
discounted by the risk free interest rate r{t) = f{t, t). In other words 

P(t,T) = %i,][e-J''*^W'^*P(t.,T)] (104) 

where, as before, [Y] denotes the average value of Y over all the stochas- 
tic variables in the time interval (t,t*). 

In terms of the Feynman path integral, eq( |104| ) yields (for measure p) 

P{t,T) = J D/p[/]e^J'''* e^[^lp(t,,r) (105) 

There are two domains involved involved in the path integral given above 
in eq.( p.05| ), namely the domain for the Treasury bonds that is nested inside 
the domain of the forward rates. These domains are shown in Figure |^. 

Although written in an integral form, the condition eq. ( P-05| ) is clearly 
a differential condition since it holds for any value of t*. Hence we take 
= t + e. The reason that we need to consider only an infinitesimal change 
for the forward rates is due to the time dependent nature of the state space 



Vt- For an infinitesimal evolution in time, the funtional integral in eq.(105) 
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Figure 7: Domains for No Arbitrage based on Treasury Bonds 

collapses to an integration over the final time variables /t+e on time slice t + e, 
that is 

P(t, T) = y" e-^(*'*)e^/^[^'^]p[/; t + e, T] (106) 

We re-write above equation in the language of state vectors, namely that 

< /i|P(t,T) >= y < /t|e-^(*'*)e-«|/i+, X /i+,|P(t + e,T) >(107) 

We have, from the completeness equation given in eq. P9|) , that 

Jt+, = y X ]t^A (108) 

Hence we have from eq.( p.O?l ) that 

</t|P(t,T)> = < /t|e-^^(*'*)e-^^(*)|P(t + e,T) > (109) 
^|P(t,T)> = e-^^(*'*)e-^^W|P(t + e,r) > (110) 

It can be verified, using the explicit representation of the zero coupon bond 
given in eq. (|6TD , that 

e+^^(*'*) |P(t, T) >= \P{t + e, T) > (111) 
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Hence, we have 

|P(t + e,r)> = e-^^W|P(t + e,r) > (112) 
^n{t)\P{t + e,T) > = (113) 

Since there was nothing special about the bond that we considered, we ar- 
rive at the differential formulation of no arbitrage, namely that any zero 
coupon bond - and consequently any coupon bond - is annihilated by the 
Hamiltonian Ti. That is 

n{t)\P{t,T) >= for all t, T (114) 

Note the similarity of above equation with the case of a single security given 
in eq. (|103| ). The role of the discounting factor is, however, very different 
in the two cases. The spot rate r is a constant for the case of a security, 
whereas in the case of the forward rates - given the difference in the domain 
for the state space at two different instances - the discounting by the spot 
rate is precisely the factor required to transform the later time Treasury bond 
P{t + e, T) to the one at an earlier time, namely P{t, T). 

7 No Arbitrage Condition for Stochastic Volatil- 
ity 

Armed with the Hamiltonian for forward rates with stochastic volatility given 
in eqs.(|77D and (|9l|) , we apply the no arbitrage conditon obtained in eq .( 11141) , 
namely 

n{t)\P{t,T) >= (115) 

or more explicity 

< ft\nit)\P{t,T) >= 7i:(t)e-^'^'^"^(*'") = (116) 

7.1 No Arbitrage for Volatililty a Function of Forward 
Rates 

Recall the zero coupon bond is given by 

P(t, T) = exp ( - /o ^ dxe'^(*'^')) (117) 
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and which yields 



S p^^^y^ _ / -/oe^(*'-')P(t,T), t<x<T 



6(l){t,x) ' ' ' [0, x>T 



118) 



We have from eqs. ( |116| ) and (^7|) that 



dxdxVoe'^'^(*'")/(^, x)D{x, x'- t)aoe'^'^(*'"')/(t, x') 

+ dxa{t,x)f{t,x) P{t,T) = 

^ a(f,a;) = ^e"<^(*'") T x'; t)aoe"<^(*'"')/(t, x') (119) 

Jo Jt 

: Condition of No Arbitrage 

Note that the no arbitrage condition given above is not contained in the 
HJM-class of solutions for the drift velocity which are all quadratic in the 
volatihty fields |jl3[; the appearance of the forward rates f{t,x) directly in 
the drift velocity emerges naturally in the field theoretic formulation, and is a 
reflection of the kinetic term in the Lagrangian for the case of / G [—oo, +oo] 

- namely - being replaced by {-^Y for / G [0, +oo]. 

We can in fact prove a more general result for the action for the case 
of when stochastic volatility is a function of the forward rates. We write the 
most general Lagrangian as 

Cgener.l = C[(j)] + J U{t,x)^ + J W{t,x) (120) 

U,W : arbitrary local functions of /(t, x) (121) 
The no aribitrage condition then yields that 

U{t,x) = W{t,x) = (122) 
In particular, a string tension term in the Lagrangian has the form 

W{t,x)cx{^y (123) 
and is forbidden by the condition of no arbitrage. 
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7.2 No Arbitrage for Volatililty as an Independent Quan- 
tum Field 



From the Hamiltonian given in eq.(^) we see that, as in the case above, 

— yields zero in the eq. ( |116|) since the zero coupon bond does not depend 
exphcitly on the volatihty field. Using the fact that 



-It dxf{t,x) 



We hence have from eqs.( |116| ), ( pT] ) and ( |124| ) 



(-1)™ e-^t dxf{t,x)^ t<x <T 
0, x>T 



(124) 



dxdx'a{t^x)G{x^x']t)a{t^x') + I dxa{t,x) P{t,T) = 



1 

~ 2 _ 

=^ a{t, x) = a{t, x) J dx'G{x, x']t)a{t, x') (125) 
: Condition of No Arbitrage 

Since there is no instrument in the financial markets at present that trades 
in volatility of the forward rates, we cannot apply the condition of martingale 
to the volatilty field, and in particular, we cannot obtain the drift velocity 
of the volatility field, namely l3{t, x), to be a function of the other fields and 
parameters of the theory. For this reason /3 has to be determined empirically 
from the market. To obtain the limit of volatility being determistic, we need 
to take the limit of p and k — 0. We then have 








(126) 


r+ - 




(127) 


r_ - 





(128) 


G{x,x';t) - 


D{x,x';t) 


(129) 



with propagator D{x,x'; t) given by eq. (ffTD . 

By a remarkable set of identities, it can be shown that the propagator 
D{x,x';t) given in eq.([71|) above is exactly equal to the one obtained in 
ref [|lOl using path integral techniques. Hence the no arbitrage condition 
obtained for a for the case of deterministic volatility using the Hamiltonian 
condition is the same as the one obtained earlier path integration. 
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Incorporating the expression for a{t,x) given in eg. ( P-25| ) into the La- 
grangian yields the final result. For notational convenience define the follow- 
ing non-local function of the volatility field 

dx'G{x,x';t)a{t,x') (130) 



The Lagrangian given by the equation above is a complete description of the 
theory of forward rates with stochastic volatility. All the parameters in the 
theory, namely the function f]{t, x) and the parameters /i, k, ^ and p need to 
be determined from market data. Due to the presence of the field v{t, x) we 
find that the Lagrangian is non-local, with the function G{x, x'; t) containing 
all the information regarding the absence of arbitrage. 

There are two further generalizations that we can make of the Lagrangian 
obtained in eq.( p.31| ), namely that the forward rate can be made positive, that 
is, / > 0, and that the propagator G{x, x'; t) can include more complex effects 
arising from a dependence on maturity of the rigidity parameter /x. 



We hence obtain 



L(t,X) = ; r- a —— 

^ ' ^ 2(1 - p2) I dt 



1 f d f _,df 

a 



2p? \dx \ dt 



8 Conclusions 

We made a generalization of the field theory model for the forward rates to 
account for stochastic volatility by treating volatility either as a function of 
the forward rates or as an independent quantum field. In both cases, the 
Feynman path integral could be naturally extended to account for stochastic 
volatihty. 



For the case of deterministic volatility, it was found in |T^ that in ef- 
fect the two dimensional quantum field theory reduced to a one- dimensional 
problem due to the specific nature of the Lagrangian. However, on treating 
volatility as a quantum field, the theory is now irreducibly two-dimensional, 
and displays all the features of a quantum field theory. 
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To exactly solve for the no arbitrage condition for stochastic volatility, 
we had to recast the condition of no arbitrage as a condition involving the 
Hamiltonian of the system. To obtain the Hamiltonian of the forward rates, 
we were in turn led to an analysis of the underlying state space of the system, 
which turned out to be non-trivial due to the parallelogram domain on which 
the forward rates are defined. The Hamiltonian for the forward rates is an 
independent formulation of the theory of the forward rates, and can lead to 
new insights on the behaviour of the forward rates. 

The model for the forward rates with stochastic volatility has a number of 
free parameters that can only be determined by studying the market. Hence 
on needs to be numerically analyze the model so as to calibrate it, and to test 
its ability to explain the market's behaviour. The first step in this direction 
has been taken in and these calculations are now being extended to the 
case of stochastic volatility. 
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Appendix 

We explicitly evaluate the propagator D{x,x';t). Note that the normalized 
eigenfunctions on the interval [t,t + Tp/j] that satisfy the Neumann condition 
of vanishing derivatives aX x = t and x = t + Tpn are given by 

^^(x) = -=cos{^^^ ^};m = 0,l,2,3, oo (132) 

yJ-FR J-FR 

which satisfy the eigenvalue equation 

19^ _ X , , X V, nin 



Hence we have 

Dixy;t) = ^^o^(x) + -^f:M^^ (134) 
2Tfr Tfr^^ (^) +1 
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and which yields, after some simphfications 



D{x, x'; t) 



2 sinh (yuTpij^ 



■[cosh(/zTp^ — — x'l) 



+ cosh (/iTp/j — //(x + x' — 2t))] (135) 
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